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Abstract 

We use the Hamilton-Jacobi theory to study the nonUnear evolutions of inhomogeneous space- 
times during inflation in generalized gravity. We find the exact solutions to the lowest order 
Hamilton-Jacobi equation for special scalar potentials and introduce an approximation method 
for general potentials. The conserved quantity invariant under a change of timelike hypersurfaces 
proves useful in dealing with gravitational perturbations. In the long-wavelength approximation, 
we find a conserved quantity related to the new canonical variable that makes the Hamiltonian 
density vanish, and calculate the non-Gaussianity in generalized gravity. The slow-roll inflation 
models with a single scalar field in generalized gravity predict too small non-Gaussianity to be 
detected by future CMB experiments. 
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I. INTRODUCTION 



Inflation scenario is a successful model to solve the problems of the standard Big Bang 
theory and explains remarkably the observational data. Quantum fluctuations of a scalar 
field are adiabatic and Gaussian during the inflation period and provide a seed for den- 
sity perturbations. The amplitudes of the perturbations freeze out when the perturbations 
stretch out to the superhorizon scale by an accelerated expansion. Inflation also gives the 
scale invariant spectrum (n^ = 1) when the perturbation modes cross the horizon. Linear 
perturbation theory is enough to explain these gravitational perturbations and temperature 
anisotropy in the early universe. 

n 

Recent WMAP observations JL] try to find a signal of the non-Gaussianity in the tem- 
perature anisotropy. The non-Gaussian signal in the CMB anisotropy might be generated 
either from the non- vacuum initial state jjor from nonlinear gravitational perturbation ^ . 
Gaussian statistical properties are completely specified by the two-point correlation func- 
tion. However, the two-point correlation function is not sufficient to describe the statistical 
properties of the non-Gaussianity, so it is necessary to investigate higher order correlations 
such as the three-point correlation for nonlinearity of a perturbation field or the four-point 
correlation for a non-vacuum initial state. Second order perturbation theory has been used 
to explain the non-Gaussianity in the temperature anisotropy In addition to second 
order perturbation theory, the Hamiltonian formalism turned out to useful to deal with 
nonlinear evolutions in the early universe and was applied to canonical quantum gravity 41 

n n 

or semiclassical gravity p, |6|| for a long time. 

Salopek and Bond in Ref. employed the Hamiltonian formalism to study the nonlinear 
evolutions of gravitational perturbations. It was also applied to Brans-Dicke theory ^ 
and low energy effective string theory [9]. Especially, the Hamilton- Jacobi theory provides 
a powerful tool to get solutions of nonlinear evolutions in the early universe through a 
generating functional which satisfies the momentum constraint equation. Even though it 
is difficult to get exact solutions of the Hamilton- Jacobi equation for general potentials, 
the large scale perturbation, which contributes mainly to the large scale structure in the 
present Universe, could be treated appropriately using the long wavelength approximation. 
The long wavelength approximation assumes that the length scale of the spatial variation 
is much longer than the Hubble radius, so it is a reasonable assumption to deal with super- 
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horizon scale perturbations. It is also known that the gauge invariant conserved quantity 



10 



exists in nonlinear perturbation theory for a superhorizon scale 

Brans-Dicke type gravity naturally emerges from the fundamental theory of particle 
physics such as string or M-theory. Although it is not clear how scalar fields couple to 
gravity, it is necessary to investigate the perturbations in the alternative gravity theory such 
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as f{(f))R type gravity as well as in Einstein gravity. Recent supernovae observation 
and WMAP results 13j imply that our universe today is in an accelerated expansion phase 
and dominated by the dark energy which has the equation of state, p/p < —1/3. It is also 
needed to consider nonlinear evolutions of such a matter component to see whether their 
existence affects the temperature anisotropy. 

In this paper, the Hamilton- Jacobi formalism will be used to study the nonlinear evolu- 
tions of inhomogeneous spacetimes in generalized gravity theory during the inflation period. 
The canonical variables will be transformed to new ones that make the new Hamiltonian 
density vanish and these new variables are constant in time for flxed spacelike hypersurfaces. 
The gauge invariant quantity, (, which is conserved in the large scale limit, is derived from 
one of new canonical variables. By introducing the nonlinear parameter, Jnl, the conserved 
quantity may be decomposed into a linear Gaussian part, C^, and a nonlinear part [l^ : 

c(x) = a(x) + ^/^l[c!(x) - (cl(x))]. (1.1) 

Non-Gaussianity is parameterized through f^i which will be constrained by observations. It 
is generally expected to be difficult to detect a non-Gaussian signal in CMB experiments for 
a single field infiation model in Einstein gravity ^3]. So the detection of the non-Gaussianity 
can constrain different inflation models. 

This paper is organized as follows. In Sec. II, we derive the Hamilton and momentum 
constraint equations in generalized gravity. The Hamilton- Jacobi equation will be obtained 
through a canonical transformation. Assuming a generating functional which satisfies the 
momentum constraint equation, we get the conserved quantity that is invariant under a 
change of timelike hypersurfaces in the long-wavelength approximation. In Sec. HI, non- 
Gaussianity will be computed using the generalized curvature perturbation on comoving 
hypersurfaces. Finally, we discuss the physical implications of the non-Gaussianity in gen- 
eralized gravity in Sec. IV. 
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II. HAMILTON-JACOBI FORMALISM IN GENERALIZED GRAVITY 



A. Hamilton equations 



The generalized gravity action to be studied in this paper is given by 



(2.1) 



where f{(l)),uj{(f)) and V{(f)) are functions of a scalar field 0. We shall confine our attention to 
the slow-roll inflation models that are described by this action. Einstein gravity is recovered 
when /(0) = I/SttG and uj{4>) = 1. Further, Brans-Dicke theory is prescribed by /(0) = 
(p/SiT, uj{(j)) = uj/8TT(j) and V = 0; the non-minimally coupled scalar field theory corresponds 
to /(0) = (1 — 87iG^(f)'^)/87TG and uj{(f)) = 1; the low energy effective string theory is given 
by /(0) = e"*^ and uj{(j)) = e~'^. We consider the Arnowitt-Deser-Misner (ADM) metric 



ds^ = {-N^ + NiN')dt^ + 2Nidtdx' + -fijdx'dx^ , 



(2.2) 



where and are a lapse function and a shift vector, respectively, and jij is a 3-spatial 
metric. The 4-dimensional Ricci scalar, R, can be written in terms of the 3-dimensional 



Ricci scalar, ^R, and the extrinsic curvature, Kij, as 



R 



^R + KiiK'^ -K^ - 



(2.3) 



The terms in the square bracket in the above formula are total derivatives or surface terms 
so that they can be integrated out in Einstein gravity. They cannot, however, be neglected 
in generalized gravity. The extrinsic curvature tensor and trace are given by 



K = Y'Ki, = Kl 



(2.4) 



2N\'"" ' dt J 

where the vertical bar is a covariant derivative with respect to 7,^. The is a generaliza- 
tion of the Hubble parameter. By varying action with respect to and 0, the momenta 
conjugate to 7,^ and are obtained as 
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From these relations, K can be written as 

1 



K 



;i + 3fi)7 
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(2.5) 
(2.6) 

(2.7) 



where vr''' = Variations of the action with respect to and lead to the Hamiltonian 

and momentum constraint equations, respectively, 
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1 



where 
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Then one finds the action of the form 



J = / d^x[n'^^ij + TTf(f) -NH- N'Hi] 



(2.8) 
(2.9) 

(2.10) 

(2.11) 



Here, and N"^ are considered as Lagrange multipliers. The evolution equations for jij and 
can be obtained by varying the action ()2.1H) with respect to vr*-' and n'^: 
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2Vt 
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-7,,7^^ (2.12) 



— vr^ — vr ' 
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B. Hamilton-Jacobi equation 



To solve the Hamiltonian and the momentum constraint equations, ()2.8|) and ()2.9|) . we 
use the Hamilton-Jacobi theory. Through an appropriate canonical transformation from 
7jj, 0, TT*-' and 'k'^ to new ones 7jj, 0, vr*-' and vf''^, we can construct a vanishing Hamiltonian 



where 



(2.14) 



(2.15) 



and the generating functional S" is a function of 0, 7jj , and 7jj. Then the canonical trans- 
formation gives the following relations 
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Finally, we get the Hamilton- Jacobi equation from the Hamiltonian constraint 
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and the momentum constraint equation 
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(2.17) 
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The momentum constraint equation implies that the generating functional S{(j),'-fij,(j),jij) 
is invariant under spatial coordinate transformations. In general, the momentum constraint 
equation does not vanish through canonical transformations as long as A^* does not vanish. 
On the contrary, the Hamiltonian constraint vanishes strongly. It is difficult to solve the 
Hamilton- Jacobi equation, (|2.17p . in general, except for special cases such as an exponential 
potential in Einstein gravity ^7j. 



C. long- wavelength approximation 

To deal with the large scale gravitational perturbations, it is reasonable to use the ap- 
proximation that temporal variations of fields are much greater than spatial variations. In 
infiation scenario, the inhomogeneous field, whose physical wavelength is much larger than 
the horizon size at the end of the infiation period, mostly contributes to formation of the 
large scale structure in the present universe and large angle CMB anisotropy. The long 
wavelength approximation assumes that the characteristic scale. A, of spatial variations is 
much longer than the Hubble radius, H~^, 0,0]: 

-dajk < ijk ^ Xph = aA > H-\ (2.19) 

where if is a Hubble parameter, and Xph and A are a physical and a comoving wavelength, 
respectively. The generating functional can be expanded in a series of spatial gradient terms 

S = + ^(2) + ^(4) + . . . _ (2.20) 

The lowest order Hamilton- Jacobi equation neglects the terms containing spatial gradients. 
In this paper we only consider the lowest order Hamilton- Jacobi equation which is sufficient 
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for dealing with the nonhnear evolution of the inhomogeneous gravitational fields. Then the 
lowest order Hamilton- Jacobi equation is 
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1 + 31]7^/V,0^""57^, 50 
We assume that the lowest order generating functional takes the following form 



+ 71/2^ _ 0. 



(2.21) 



(2.22) 

where ^ is a new canonical variable involving 7jj and 0, and /3 is a constant which carries a 
dimension. For Einstein gravity, W can be interpreted as a locally defined Hubble parameter 
if we taUe , ^ -^O \k. Th. ,ene.a«„, func«„„a. au.o,na«ca„y sa«s«es the .o.e.u. 
constraint equations if A^* = 0, and this will be discussed in the next section. 

It is convenient to factor the 3-spatial metric 7jj into a conformal factor and a conformal 
3-spatial metric hij with the unit determinant det(/ijj) = 1: 



7ij(t,x) =7^/3(t,x)/i,j-(x). 



(2.23) 



The gravitational waves are related to the hij. Then, from Eq. ()2.1(j|l . the conjugate momenta 
for 7jj and are 
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(2.24) 
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(2.26) 



which is traceless, 'jijdW/d'-fij = 0. If we decompose vr*-' into the trace vr'^ and the traceless 
part vf*-' , they are given by 
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(2.27) 



As long as we do not concern about the gravitational radiations, W is assumed to be 
independent of hij, so vf*-' can be set to zero. Then the Hamilton- Jacobi equation becomes 
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and the evolution equations for 7 and are 
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(2.29) 
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D. Solutions of Hamilton- Jacobi equation 

Although the Hamilton- Jacobi equation (j2.28|) is difficult to be exactly solved for general 
potentials, it can be solved for some special cases in generalized gravity. 

1. V = 3^2/2 

If the potential is given by V{(f)) = 3/^2/2(0)^ then the Hamilton- Jacobi equation takes 
the form 

2 



dW 



where 



2/ • 



As ^2 > 0, Eq. ()2.3H) can be exactly integrated to yield [17 1 



W{(j)) = cosh[± J a{(f))d(j) + C], 



(2.31) 



(2.32) 



(2.33) 



where C is an integration constant. The solution of the form ()2.33|) is available for a constant 
potential in Einstein gravity or a massive scalar field potential for the generalized gravity 
with /(0) oc 0. 

2. V = or general potentials 



With the identification of W = q and = t, the Hamilton- Jacobi equation may be 
interpreted as a time- dependent inverted oscillator with a unit mass, frequency a, and 



energy a p: 



ifowV _ 1 

2\d^) "2 
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FIG. 1: The inverted oscillator for p = 0: g 




= 0. 



where = V/ {SP"^ p). In fact, this inverted oscillator has a 0-dependent energy and 



First, in the case of p = corresponding to the Brans-Dicke gravity or low energy effective 
string theory, the solution can easily be obtained 



where Wo(0i) is an initial value at 0,. When a is non-square integrable, the solution ()2.35|) 
for the upper (+) sign grows to ±oo as goes to oo depending on the sign of Wo{(l>i). 
This corresponds to the downward motions in Fig. 1. Whereas, for the lower (— ) sign, 
the solution approaches an attractor regardless of Wo{(f)i) as goes to oo. The solution 
approaches the attractor regardless of Wo(0j) for the upper sign, but it diverges to ±cxo 
for the lower sign depending on the sign of Wo(0j) as goes to — oo. On the other hand, 
for a square integrable a, as goes to ±oo, the solution ()2.35|) approaches finite values, 
Wo(0i) exp[± /^°° ad(f)], not necessarily attractors, which correspond to the upward motions 



For general potentials with p 7^ 0, we have a similar picture as shown in Fig. 2. Not 
only the energy but also the curvature of the potential depend on 0. In the analogy of an 
oscillator, as shown in Fig. 2, the total energy line moves up or down depending on a^p 
and the 0-dependent curvature narrows or widens the parabola. As for the p = case, 
W either approaches to attractors or diverges to ±00, depending on the sign of Wo{(f)i) and 
the behavior of a and p, as |0| grows. This behavior of W for general potentials can be 



curvature (spring constant) of potential. The W = ±y^ are two fixed or stationary points. 




(2.35) 



in Fig. 1. 
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FIG. 2: The inverted oscillator for p / 0: i - ^a^W^ = -\a^p. 

calculated numerically [isj ]. 

To find approximately an analytical solution for general potentials, we can write W{(j)) 
as W{(t)) = Wo{(f))Z{(f)), where Wo^cp) is a homogeneous solution ()2.35p . Then the Hamilton- 
Jacobi equation reduces to the equation for Z{(f)): 

where we have used Eq. (|2.34|) for p = 0. Equation (j2.36j) is a nonlinear equation for Z{(j)). 
Assuming that Z{(j)) is a slowly varying function of 0, we can introduce a small parameter 
6 to indicate smallness of the nonlinear terms and rewrite Eq. fl2.36p as 

JlnWo ^dZ (dZy a^p 

The parameter 6 will be set to one in the final result. Now we expand Z{(j)) in a series of S 

Z{(j)) = + + (52^(2) + . . . . (2.38) 

Substituting Eq. ()2.38j) into Eq. ()2.37j) and comparing terms of the same powers of S, we 
obtain the following equations up to first order 



5°: 

d\nWo ^(Q^dZ^^'^ a^p 
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d(P \ d(t) d(t) ) \ del) ) ^ ' 



The solutions at each order can easily be found 

,2^ \V2 



a p 



' W^d\nWo/df 

^(1) ^ \_ f {dZ^'Vdct>) 

2Z(o) J dW^/dcj) 



2 

-ci0. (2.41) 



On the other hand, for a rapidly varying function 2'(0), we can repeat the same analysis 
and write Eq. ()2.36|) as 

2 ..r 9, 



With the series expansion ()2.38p for Z{(j)), we obtain the following equations up to first order 

^ + ^Z<«..0 (2.44) 

(70 (70 

These equations can be simply integrated to give 



Z« = -/^Z(°)rf0. (2.45) 



The validity of the approximation in Eqs. p.37|) and p.42|l for specific models will be 
discussed in a future work Q]- 



E. conserved quantities 

In Eq. fl2.22|l . W depends on 0, 7jj, and 6', where 6' is a new canonical variable that has 
the conjugate momentum 

vr^ = -2/?y/V^/^(0)^. (2.46) 
11 



Then the momentum constraint equation ()2.9p reduces to j3] 

Hi = TT^diO. (2.47) 

Since the new variable is chosen to make the new Hamiltonian density vanish, the new 
Hamiltonian contains only a contribution from the momentum constraint 

H = J d^xN'TT^diO. (2.48) 

The Hamilton equations for the new variable are 

e - N'diO = 0, TT^ - diiN'-n^) = 0. (2.49) 

In general, the new canonical variable needs not to be constant in time, but if the spacelike 
hypersurfaces are chosen such that the shift vector, iV\ vanishes, then 6 and tt^ are constants 
for fixed spatial coordinates Q|. 

It is known that a gauge invariar it q uantity can be derived by taking the spatial gradient 
of an inhomogeneous quantity llQ^ . By taking the spatial gradient of the logarithm of 
the canonical new variable ()2.46p . we can obtain the gauge invariant quantity 



3d\nf 5 , dW 
+ TT-ln 



(2.50) 



2 d(f) d(f) de 

To calculate the last term in the above equation, we differentiated the Hamilton- Jacobi 
equation ()2.28|) with respect to 6: 

With this relation, Eq. ()2.5()|1 now becomes 

di\n\Tx'\ = ]-di\n-i + ^d4, (2.52) 
2 (p 

where we have used Eqs. fl2.29p and ()2.30|) . This quantity is similar to the gauge invariant 
quantity C = \I' — {3H/(p)6(f) in the linear perturbation theory, which is conserved at the 
superhorizon scale. Here \1/ is a Newtonian gravitational potential. We define the generalized 
gauge invariant quantity Q in nonlinear theory, which is conserved in the large scale limit 

Q = ldi\n\TT'\ = \dM^ + ^d4. (2.53) 
3 6 30 

Although Ci is defined on the uniform energy density hypersurfaces in general, it coincides 
with a curvature perturbation T^j, which is defined on the comoving hypersurfaces in a scalar 
field dominated universe and also conserved in the large scale limit. Q and TZi are proved to 



be conserved in the large scale limit in Refs. 
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III. NON-GAUSSIANITY IN GENERALIZED GRAVITY 



Non-Gaussianity in CMB might be generated by a non-vacuum initial state |yj or a 
nonlinear perturbation P], even though the initial perturbation is Gaussian. Although the 
non- vacuum initial state gives the zero three-point correlation function, the relation between 
the four-point and the two-point correlation functions, which obey the Gaussian statistics, 



\ T 



((^)) = 3((f)V, 



(3.1) 



is no longer satisfied P|. Whereas the nonlinear gravitational perturbation leads to a non- 
zero three-point correlation function. In this paper we shall focus on the non-Gaussian signal 
in CMB only from nonlinear perturbations. To show the non-Gaussianity by the nonlinear 
perturbation, the gravitational potential $ may be decomposed into a linear part and a 
nonlinear part with a nonlinear parameter /atl |1J] 



$(x) = $^(x) + /^,[$i(x)-(<l.i(x))]. 



(3.2) 



Here, $l(x) is a linear Gaussian perturbation that has the zero expectation value = 0. 

With this definition, the non- vanishing component of the $(k)-bispectrum, which is the 
Fourier transform of the three-point correlation function in the coordinate space, is Q] 



(<I>Z.(ki)<l>L(k2)<l'7VL(k3)) = 2(27r)35(3)(k^ + k2 + k3)/^z.P*(ki)P$(k2), 



where P$(k) is the linear power spectrum given by 



(3.3) 



(<l>z.(ki)$i(k2)) = (27r)3p$(ki)5(=^)(ki + k2). 



(3.4) 



and 



$JVL(k) = /. 



NL 
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.$,(k + k')$I(k')-(2vr)35(=^)(k)(<l>i(x)) 



(3.5) 



It is known that /^vj;^ should be larger than order unity to be detectable by CMB experi- 
ment. But the single field inflationary rnodel gives too much small value of f^i oc 0{e,T]) 
where e and rj are slow-roll parameters ^. Thus, if the non-Gaussianity is detected, it can 
constrain inflationary models. However, it would be interesting to calculate non-Gaussianity 
in generalized gravity theories. 
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We follow the method in Ref. [3] to calculate the nonlinear curvature perturbations on 
the comoving hypersurfaces. The 3-spatial metric jij can be written as 



7ij(t,x) = a^{t,x.)hij{x.) 



(3.6) 



where a(t,x) is a local expansion factor, and the conformal 3-metric, hij{x.), is independent 
of time and has the unit determinant, det(/ijj) = 1. Then the local Hubble parameter H 
takes the form 



H = 



1 a 
Na 



3 



(3.7) 



If we choose as a time coordinate on the comoving hypersurfaces, we can obtain lna(0, x) 
from Eq. (O:^ 



In a{(f>, x) — In a{(j)o, x) 



<f>0 



d(f)NH. 



(3.8) 



We use the notation (50)7^(lnao, x) = 9j0(lnao,x) = 0(lnao,x) — 0(lnao,Xo) for the 
scalar field fluctuation on the spatially flat hypersurfaces {dilna = 0) and TZs^{(f)o,x.) = 
dilna{(j)Q,x.) = lna(0O)X) — lna(0O)Xo) for the curvature perturbation on the comoving hy- 
persurfaces {di(f) = 0). The scalar field fluctuations on the spatially flat hypersurfaces are 
needed to transform to the curvature perturbation on the comoving hypersurface. Then the 
curvature perturbation on comoving hypersurfaces is given by 

r4>o+iS4>)n 



T^5(j, — — 

Using Eqs. (IT^ and TH^ . we get 7^^ 

r4>o+{S4i)Tz 



</>o 



d(pNH. 



(3.9) 
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-1 
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where 



D(0) 



cu(l + 3fi) 

7 ■ 



(3.11) 



By expanding the nonlinear relation between 7^^^ and {dcfyn up to second order 
obtain a nonlinear curvature perturbation, IZ^^ = IZl + 'R-nl, where 



we 
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1 
G 
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dDdlnW B'^lnW (dwVd'^lnf 
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+ 



(3.12) 
(3.13) 
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where 



= D + 



d\nfd\nW' 



(3.14) 



Finally, from Eq. ()3.2j) . we obtain the nonlinear parameter 



NL 



5G 



dDdlnW d^lnW fdW\ d'^lnf 



(3.15) 



where we have used the relation between $ and 7^ as $ = ^TZ during the matter dominated 



era. If we define the slow-roll parameters in generalized gravity theories by 

H 6 I dH 



N\NJ Hcj) 
using again Eqs. (j2.29|) and (j2.3U|) . we obtain 

1 dlxiW 



(3.16) 
(3.17) 



G 



dlnf j d'^lnf ifdlnfY d\n D d\n f 



+2 



dlnf 1 d'^W 



^0 W 

1 d\nW 



(3.18) 



G 



D - 2^}^] + 2D (^] 



^f/51n/\ ^d\nfdlnD)d\nW ^^dlnf 1 d'^W 



(3.19) 



/ J W (902 

The nonlinear part of the curvature perturbation can be written in terms of the slow-roll 
parameters, e and r], as 

7^7VL = ^(r?-e)(7^L)2. (3.20) 

Hence, the nonlinear parameter, Jnl, for generalized gravity theories takes the form 

fNL = l{v-e). (3.21) 

Note that the nonlinear parameter ()3.21|) in generalized gravities has the same form as in 
Einstein gravity. The WMAP result gives a constraint on /jvl by —58 < Jnl < 134 at 95% 
CL 0. Since the slow-roll inflation implies that |e|, \ri\ <^ 1, non-Gaussianity in the single 
scalar field inflation model is difficult to be observed by CMB experiments. 
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IV. CONCLUSION 



It is expected that nonlinear perturbations may be responsible for a non-Gaussian signal 
in CMB experiments. The Hamilton- Jacobi theory can provide a useful and convenient tool 
to deal with nonlinear perturbations. In this paper we derived the Hamilton- Jacobi equation 
in generalized gravity theory. Through a canonical transformation, a new set of canonical 
variables was chosen that could make the new Hamiltonian density vanish. The conserved 
quantity, Q = dilna — {NH/(p)di(j), was obtained in the large scale limit by using the fact 
that the new canonical variable is constant in time on fixed spacelike hypersurfaces. The 
(i could be regarded as a generalization of the gauge invariant quantity ( = — {?>H/^)5(f) 
in the linear perturbation theory. In the long wavelength approximation, we found the 
exact solutions to the lowest order Hamilton- Jacobi equation for special scalar potentials 
and introduced an approximation scheme for general potentials. 

The non-Gaussianity from nonlinear density perturbations was parameterized with a 
nonlinear parameter Jmli which is an expansion parameter of the gravitational potential in 
Eq. (j3.!2|) . The nonlinear parameter can be measured by CMB observations. Any detection 
of the non-Gaussianity of CMB may put a strong constraint on inflation models. Especially, 
Jnl predicted by the single field slow- roll inflation in Einstein gravity is too small to be 
detected in CMB experiments through a non-Gaussian signal. Nevertheless, it would be 
interesting to investigate the non-Gaussianity in generalized gravity theory. We have found 
that even in generalized gravity theory, the nonlinear parameter for the slow-roll inflation 
takes the same form /tvl = (r^ — e)/2 as in Einstein gravity. Hence, the slow-roll inflation 
models based on Brans-Dicke theory, non-minimally coupled scalar field theory, and low 
energy effective string theory, et al as well as Einstein gravity will have the same order of 
non-Gaussianity for CMB. This has a physical implication that slow-roll inflation models 
with a single scalar field in such alternative theories of gravity will be ruled out if a signal for 
non-Gaussianity is observed in future CMB experiments. Therefore, the non-Gaussianity 
would require multi-field inflation models or a different generating mechanism for density 



perturbations 



2l| in generalized gravity. 
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